Nonlocal impedances and the Casimir entropy at low 

temperatures 

V. B. SvetovoJ] 

MESA+ Research Institute, University of Twente, 
P.O. 217, 7500 AE Enschede, The Netherlands 



& 



R. EsquiveJ 

Instituto de Fisica, Universidad Nacional Autonoma de Mexico, 
Apartado Postal 20-364, DF 01000 Mexico, Mexico 
(Dated: February 1, 2008) 

Abstract 

The problem with the temperature dependence of the Casimir force is investigated. Specifically, 
the entropy behavior in the low temperature limit, which caused debates in the literature, is 
analyzed. It is stressed that the behavior of the relaxation frequency in the T — > limit does not 
play a physical role since the anomalous skin effect dominates in this range. In contrast with the 
previous works, where the approximate Leontovich impedance was used for analysis of nonlocal 
effects, we give description of the problem in terms of exact nonlocal impedances. It is found 
that the Casimir entropy is going to zero at T — > only in the case when s polarization does not 
contribute to the classical part of the Casimir force. However, the entropy approaching zero from 
the negative side that, in our opinion, cannot be considered as thermodynamically satisfactory. 
The resolution of the negative entropy problem proposed in the literature is analyzed and it is 
shown that it cannot be considered as complete. The crisis with the thermal Casimir effect is 
stressed. 
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I. INTRODUCTION 

An attractive force between uncharged metallic plates, predicted in 1948 by Casimir Q], 
is one of the most striking macroscopic manifestations of quantum vacuum. Recently this 



force became a subject of systematic experimental investigation 
force between ideal metals at zero temperature if, 



The 



K ' 2llC (1) 



240a 4 ' 

depends only on the separation a and fundamental constants. In reality the force is measured 
at finite temperature between deposited metallic films, which have finite conductivity and 
roughness. Correction to Eq. due to finite conductivity can be as large as 50% for small 
separations a ~ 100 nm. Contribution of the finite temperature to this correction is not 
large but caused a lot of controversy in the literature (see Ref. Q| for a recent review). The 
essence of the problem lies in the classic contribution to the Casimir force, which dominates 
at large distances between plates or at high temperature. Calculations made for ideal metals 
at finite temperature jlll . \v\ showed that s- and p-polarized modes of electromagnetic field 
gave equa 



contributions to the force. At the same time the Lifshitz theory of fluctuating 
fields Q,Q predicted zero contribution for s-polarization. For the first time the problem 
was recognized many years ago. For reconciliation of the results Schwinger, DeRaad, and 
Milton (SDM) [15] proposed a special prescription to be used with the Lifshitz formula: 
one has to take first the limit e —>■ oo for the metal permittivity and only then allow the 
frequency uo go to zero. Modern calculations concerned with nonideal metals were confronted 
with the problem again. 

Different approaches to resolve the problem have been proposed in the literature, which 
resulted in different temperature corrections to the Casimir force. Bostrom and Sernelius 
used the Lifshitz formula with the Drude dielectric function and found that s-polarization 
did not contribute in the classical limit (n = term in the Lifshitz formula) independently 
on the Drude parameters. In this approach there is no continuous transition to the ideal 
metal case and the predicted temperature correction is in contradiction with the Lamoreaux 
experiment ;2|. However, physically this approach is well motivated since the Drude dielectric 



function is working especially well at low frequencies. Bordag et al. 



171 ] used the plasma 



model dielectric function, for which e at low frequencies increases faster ( uj 2 ) than for the 
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Drude model They found that s-polarization gives finite contribution in the classical 

limit, which coincides with the ideal metal result when the plasma frequency uj p is going to 
infinity. The temperature correction in this approach is very close to that for the ideal metal 
and negligible at small separations between plates. A weak point of this approach is that 
no known material behaves at low frequency according to the plasma model. Svetovoy and 
Lokhanin proposed to use SDM prescription for the n = term in the Lifshitz formula 
for real metals also. Later it was shown that this prescription follows from very general 
dimensional analysis of the classical contribution to the force if one demands continuous 
transition to the ideal metal case. The temperature correction happened to be small but 
observable at small separations between bodies. 

A new round of discussion has started when a thermodynamical problem connected with 
the Casimir free energy has been revealed \'2(\ . The idea was to use the Nernst heat theorem 
as a guiding principle to choose between different approaches to the temperature correction. 
According to this theorem the entropy has to go to zero in the limit of zero temperature. 
It was noted j^j] that the Drude relaxation frequency u T vanishes with T and, therefore, 
the plasma dielectric function is realized at T — > 0. In this case the leading term in the 
temperature correction is ~ T 3 and the entropy is safely going to zero as S ~ T 2 . Two 
other approaches predict the leading term in the correction ~ T and finite entropy at T = 0, 
positive and negative for the approaches H and Q, respectively. However, the following 
analysis revealed that the situation is not as simple. The anomalous skin effect was shown to 
be important for the temperature correction at low temperatures j^jj. With the use of the 
Leontovich impedance for the anomalous skin effect it was demonstrated that the entropy is 
going to zero only if SDM prescription is used for the n = term. On the other hand, it was 



noted 



22j and expressed later more clearly |23J that any real material contains a number 



of defects, which are responsible for the residual resistance at T = 0. Equivalently it means 
that uj T becomes very small but finite at T = 0. It was shown, that in this case, the entropy 



disappears at sufficiently low temperature 



Therefore, again we have a confusing 



situation where each approach has its own reasoning. 

We would like to emphasize that at low temperatures the anomalous skin effect plays an 
important role and should be taken into account in any reasonable calculations. Because u T 
decreases fast with the temperature, at sufficiently low temperature inevitably the mean free 
path I = vf/uj t (T) for electrons becomes much larger then the field penetration depth S. 
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When this happens the relaxation frequency does not play a physical role any more. Instead 
of uj t the physical significance gets the other frequency, Q = (vp/c) u p , which is often used 
as a characteristic frequency of anomalous skin effect. For this reason the question, does u T 
go to zero or have some residual value at T — > 0, becomes unimportant. This ideology was 
developed in Ref. [21] in context of the thermal correction to the Casimir force. For the 
description of the anomalous skin effect the Leontovich impedance was used there. 

The approximate Leontovich impedance was used for calculations (21], 24, 2^ (see addi- 
tional discussion in Refs. l26l. 1271 128IH . The approach similar to the Leontovich impedance 
was developed also in Refs. |2iJ, |3i|. This impedance describes well the propagating electro- 
magnetic field, but it was not clear why in the local limit it gives the result different from 



the dielectric function approach 



2l|. In Refs. [3l|, |32j, |33J] it was demonstrated that the use 



of the exact impedances is in agreement with the dielectric function approach. It became 
clear that the point of contradiction is the transverse momentum, which is neglected in the 
Leo^icn i.pedance Q Q Q. In o„ r paper Q a g ene r a, approach to the nonioca, 
impedances was developed for applications in the Casimir force calculations. It was shown 
that for real metals both contributions in the force from propagating and evanescent fields 
are important. The propagating fields can be described well by the Leontovich impedance, 
but the same is not true for the evanescent fields. The latter ones should be described by 
more general impedances, for which dependence on the transverse momentum cannot be 
neglected. Explicit expressions for these impedances were presented in Ref. [.36']. 

It is important to notice that the relevance of spatial dispersion effects depends on the 
separation between the slabs, being more important at short separations. For two Au slabs, 
at a separation of the order of the plasma wavelength of Au(130 — 140nm), the difference 
between the local and nonlocal calculation is 0.2%, and can be significant when experimental 
errors of the order of .5% are claimed |8[. For the hydrodynamic model ref. [32jand ref. 
[37! ] give the same results, using ajiielectric function valid in a wide range of frequencies not 
only at the infrared as stated in 



Inadequacy of the Leontovich impedance forced us to reconsider the result of Ref. [2jj for 
the entropy behavior in the low temperature limit. In this paper we calculate analytically 
the temperature correction to the Casimir free energy using the general approach to the 
nonlocal impedances js^J. 

The paper is organized as follows. In Sec. HI we separate the temperature dependent 
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part of the free energy and transform it to the form convenient for calculations. In Sec. 
IIHI the nonlocal impedances at low temperature are discussed. In Sec. IIVI we give analytic 
expressions for for the free energy in two limit cases. The entropy behavior at T — > and 
discussion around are given in Sec. In the last section we present our conclusions. 



II. TEMPERATURE DEPENDENT PART OF THE FREE ENERGY 

The Casimir force at nonzero temperature between plates made of real materials is given 
by the Lifshitz formula ^| . For the free energy J 7 (a, T) this formula can be presented in 
the following form 

oo 

t) = £L Yl ' I d y y{ ln I 1 - r s *) e ~ y ] + (*■-- . ( 2 ) 

n=0 { 

where £ n are the dimensionless Matsubara frequencies defined with respect to the charac- 
teristic frequency u a 

y- (n >. 1-nkT c 

?f» = — , Cn = Wa = — ■ (3) 

In Eq. (J2J) r s and r p are the reflection coefficients for s and p-polarizations, respectively. 
The integration variable y is defined via the physical values as 



y = 2a^Q/c 2 + q\ (4) 

where q is the absolute value of the wave vector along the plate. 

The problem with the thermal correction comes from the n = term in Eq. (121) . w hich 
will be denoted as JF (a,T). There is no agreement between different authors [ly, IrJ, llj 
what is the reflection coefficient r s (0, y) in this term. The n = term describes the classical 
contribution to the free energy, which dominates at large separations or high temperatures. 
Without loss of generality it can be parameterized as 

kT 

F (a,T) = -a—((3). (5) 

Here a is a dimensionless function of material parameters and separation a. This form of 
JF follows from a simple dimensional analysis [19] in the classical limit (not considering the 
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Plank constant K). Different approaches to the temperature correction problem give different 
values of a. This value will be kept arbitrary in the calculations and will be specified only 
for the discussion of the final result. 

We are interested in the temperature dependent part of the free energy, which is respon- 
sible for the entropy. To separate the temperature independent part, let us rewrite the free 
energy in the following form 

fc _ 00 

T{a, T) = Jo (a, T) + —^7^ V [G s (nr) + G p (nr)] , (6) 

n=l 

where the functions Gi (nr) (i = s,p) are defined as 

00 

Gi (nr) = J dyy In [l - rf (nr, y) e~*] (7) 

nr 

and nr was introduced instead of £„. The parameter r, 

r = 7^ — 1 kT e tf — — — fku a , (8) 
T e// ' ;/ 2a V ; 

is a dimensionless temperature. It is convenient to rewrite the sum in Eq. (JHJ) using the 
Abel- Plana formula 



1 00 
^ ^ o T ^ fGi(T + itr) 



e 2nt _ I 



dt 



^^G l (nr) = i- y G i (x)dx + ^ -G^r)- J G l (rt)dt-2lm J 

n=1 L 

'(9) 

The first term on the right hand side does not depend on temperature, but all the other 
terms describe the temperature correction. The temperature dependent part of the free 
energy AT (a, T) = T (a, T) — T (a, 0) can be presented then in the following form: 



kT 

AT(a,T)=T (a,T) + —^ 



1 00 
\G S (t)- J G s (rt)dt-2lm J G dt + (s -» p) 







(10) 

It is important to see clearly which frequencies give the main contribution to AT(a,T). 
Indeed, the most important contribution to the temperature independent part comes from 
the Matsubara frequencies ( n ~ u a or nr ~ 1. The same is not true for A J 7 (a, T). As one 
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can see from Eq. (jlOj) the important values of the dimensionless frequency £ = rt are of the 
order of 1 or Q ~ ™ a . We are analyzing the temperature behavior in the low temperature 
range, where r«l. Therefore, frequencies much smaller than the characteristic frequency 
u a give the main contribution to the temperature dependent part of the the free energy. 



III. NONLOCAL IMPEDANCES AT LOW FREQUENCIES 

As was said in Introduction, at low temperatures the importance of the anomalous skin 
effect significantly increases. Description of this effect is given within the theory of nonlocal 
interaction between the electromagnetic field and a metal. In this theory, the reflectivity of 
the metal is described by the surface impedances. The impedances are connected with the 
nonlocal dielectric functions by the general relations j3] 

oo 

Z a (u,q) = -- / - z — , (11) 

tt c J [uj l jc 1 ) e t — k 



— oo 

oo 



iu f dk z 



q 2 k 



+ 



2 



(12) 



— oo 



{uj 2 /c 2 )s 1 {uj 2 /c 2 )e t -k 2 \ ' 

where k = \J q 2 + k 2 is the wave number, e t (k, uj) and e\ (k, to) are the nonlocal dielectric 
functions describing the material response to transverse and longitudinal electric fields, re- 
spectively. These equations are true independently on the particular model used for obtain- 
ing e t and E\. The dielectric functions can be found, for example, by solving the Boltzmann 
kinetic equation 



j. The Boltzmann approximation is valid in the range uo < uj p and 
q < kp, where kp is the Fermi wave number, and is appropriate for our problem. 

The impedances Eq. tllljl and Eq.()12j) are analytic functions in the upper half of the 
complex frequency plane and can be written at imaginary frequencies u = i( using the 
analytic continuation. Explicit form of the dielectric functions along the imaginary axis 
is 

3 w-arctanv 
£liCv) = 1+ UCT^y fl{v) = V 2 ' v + MQ{v- arctan,) ' (13) 

£t(c ' w) = 1+ c7crSy' ^) = ^H+( i+ <> arcta H> ( i4 ) 



7 



where Vp is the Fermi velocity. The range of the anomalous skin effect corresponds to 
large values of v. When the Casimir force is calculated, k is restricted by the condition 
k ^ q ~ l/2a. On the other hand, the denominator in Eq. is small and the condition 
v ^> 1 will be fulfilled at sufficiently low temperature. In this limit the dielectric functions 
behave as 

One can immediately see that the relaxation frequency falls out from the dielectric func- 
tions. The longitudinal function, ei, does not depend on frequency at all, but k dependence 
describes the Thomas-Fermi screening of the longitudinal electric field. In the transverse 
function, e t , the term vpk plays a role of the relaxation frequency. The surface impedances 



corresponding to the functions (fTU|) . (JT7J) were found in Ref. 36]: 



Z.(C,q) = -F(b), (18) 
cq 

Z p ((,q) = ^= C 7 ^ + -G(b), (19) 
V3 (u p cq 

where the functions F (b) and G (b) are defined as 

F(b) = - I dx C ° Sh X G(b) = - f dx Sinh X (20) 
7r J cosh 3 x + b 3 ' vr J cosh 3 x + b 3 ' 

o o 

i= i/*r«£cr (21) 

g V 4 c ^/ 

The asymptotics for large and small values of b are 

F(6) = 1 + 0(b 3 ), G{b) = \ + 0(b 3 ), b < 1, (22) 
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m = + 0(r3) - G(6) = + 0(r3) - b>1 - (23) 



The Leontovich impedance for the strong anomalous skin effect is reproduced at finite 
frequency in the limit q — > when 6 ^> f : 

WO.^KJ^K)^^) . (24) 

However, the most important contribution to the Casimir force give finite values of q ~ 
l/2a and the limit i > 1 inevitably will be broken at some sufficiently low frequency 
(temperature). When ( is so small that b <C 1, the impedance Z s approaches the local limit 
Zs{q,C) = C/ c <?) which does not depend on u T . This is in contrast with the Leontovich 
impedance, which behaves in the local limit as Z s (() = e^Q^ 1 ^ 2 — > yj (u T /ujp. It clearly 
depends on the value of uj t . Indeed, the reason for this is the non-dependence on q of 
the Leontovich approximation. The impedance Z p also behaves very differently from the 
Leontovich impedance in the limit b <C 1, but in contrast with Z s it is significantly nonlocal. 
This is because for b <C I the main contribution in Z p gives the first term in Eq. (fTTJj) 
responsible the Thomas-Fermi screening. 

Let us discuss now the temperature range where Eqs. (JTSJl and (JT§j) are true. The main 
condition is that the parameter v in Eq. (|15|) has to be large. The minimal value of the 
wave number is k = q ~ l/2a. The important frequencies contributing to the temperature 
dependent part of the free energy Eq. fflTjj) are £ ~ 2irkT/h. We assume that the relaxation 
frequency uo T decreases with temperature faster than linearly and for this reason it can be 
neglected in the denominator of Eq. (|15|). At very low temperatures this assumption can be 
broken due to residual resistivity, but in this case v will be certainly large. Therefore, the 
value of v will be much larger than 1 if 

When this condition is met, the impedances Eq.([T%|). Eq.([T§|) can be used independently 
on the value of b. However, for large b, when the Leontovich impedance ()24j) can be used, 
the condition on the temperature is relaxed. This is because for b ^> 1, the wave number 
k = gcoshx ~ qb ^> l/2a. In this limit the condition t)>l means 
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kT«^^ (26) 
4 V vr c v ' 

This is the restriction on the temperature used in Ref. [2JJ] . Therefore, the Leontovich 
impedance can be used in the temperature range Huo a (vF/c) <C 2nkT <C a/ '3/ lQirfkOp^F / c) . 
When the temperature is going down and obeys the condition 2-akT <C hu a (vF/c), the in- 
dependence of the impedances becomes important and one has to use Eqs. fITHj) . (fT^j) . and 
Eq. d2U) instead of Eq.(j2U). 



IV. EVALUATION OF THE FREE ENERGY 

The temperature dependent part of the free energy is defined in Eq. (jlOJ) . where the 
functions G s>p are given by Eq. (|7J). The reflection coefficients be expressed via the 

impedances as 

„ _ Z s0 - Z s ^ _ Z p0 - ^ p ^ 



where Z s0 and Z p0 are the "impedances" of the plain wave defined as the ratios of the electric 
and magnetic fields in the wave: 

Z t0 =4r> Z p° = ^?> ko=V( 2 /c 2 + q*. (28) 

CKq (, 

Let us first calculate the function G s (r). Using Eq. (|TH|) for the impedance Z s the 
reflection coefficient can be written in the following form 



1 " -HHFiA/y^^) 

(29) 



\/y 2 --r 2 



where we introduced the parameter A similar to that in Ref. [21], which is defined as 

67T C W„ 



A = — —2 T ) ■ ( 30 ) 
In Eq. near the lower integration limit y ~ r the reflection coefficient tends to 

1 because the argument of the function F is in the Leontovich region (b ^> 1) where 
(y/vV - r 2 )F(A/y/y 2 - r 2 ) ~ r/A < 1. For this reason the contribution to the inte- 
gral from the region nearby the lower limit will be ~ r 2 . One can neglect this contribution 



10 



changing the lower limit by zero and neglecting r 2 in \Jy 2 — t 2 . Introducing the integration 
variable x = y/A one finds 



oo 

G s (r) =A 2 j dxx In (l - r 2 s (l/ x)e~ Ax ) + 0(r 2 ), (31) 



where 

l-F(l/x) 



(32) 



l + F{l/x) 

This integral can be analyzed in two limit cases A ^ 1 and A ^> 1. The former case is 
realized for extremely low temperatures (T <C 0.1° K for a > 100 ram) and should be used 
to check the behavior of the entropy at T — > 0. At all realistic temperatures the latter case 
is realized. Let us consider this case first. At A ^> 1 the main contribution to the integral 
Eq.pip comes from the region x ~ 1/A. Then the argument of F(l/x) is large and we are 
in the Leontovich impedance region, where F(l/x) ~ 4x/3y / 3- This situation was already 
described in Ref. 2l| and the result can be written immediately: 

G s (r) = -C(3) + C(3) + 0(1/ A 2 ), A » 1. (33) 



Here an additional factor 4/3v^3 in comparison with Eq. (29) in Ref. [21[ takes into account 
different definitions of A used in this paper. Two other terms in Eq. (fTUj) can be easy 
calculated with the help of Eq. (J22J) • 

In the opposite limit A <C 1 the situation is different. In this case, the important values 
of x in the integral (|3ip are x ~ 1 and the Leontovich approximation is no longer valid. In 
this limit the exponent can be changed by 1 and the integral in Eq. (}3T]) is just a number 
that can be found numerically substituting expression for F[l/x) from Eq. ()20|) into Eq. 
()32|) . The result will be the following: 

G S ( T ) = -0.0938A 2 + 0(A 3 ), A < 1. (34) 



The important difference of this expression from that found in Ref. 21] (see Eq. (23) 
therein) is that G s — > when A is going to zero instead of the finite value G s — > — C(3). 
The reason for this change of the behavior is the reflection coefficient. When the Leontovich 
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approximation is used in Eq.([32|) r s — > 1 at A — > but the use of the exact impedance (fTHJl 
gives r s — > 0. 

Now let us find the function G p (t). As in the case of G s (t) one can neglect r 2 in \/y 2 — r 2 
and change the lower integration limit by zero. In the low frequency (temperature) range 
the first term in Eq. (fTTj|) dominates and the reflection coefficient can be presented as 

Typically the reflection coefficient for p-polarization is approaching 1 in the low frequency 
region. Small correction in Eq. (|H5jl appears as a nonlocal effect connected with the Thomas- 
Fermi screening. With this r p the integral in Eq. (JJJ) is easily calculated and for G p {r) one 
finds 

C ,( T ) = -C(3)(l--^). ,36, 

It does not depend on r at all and holds true in both limits of large and small A. The same 
conclusion was made in Ref. j^] but without the Thomas-Fermi correction. 

Now we are able to present the final expressions for the temperature dependent part of 
the free energy in the limits A -C 1 and A ^> 1. Calculating the integrals in Eq. (|10|) using 
the functions Eq.([34p. Eq. (j36|) and the definition of the n = term in Eq.fjSJ) one finds in 
the limit of small A: 



AF(a,T) 



kT 



57TGT 



-aC(3) + -C(3) ( 1 



+ 0.0146A 2 + 0(A 3 



A < 1, (37) 



where the first term, containing a, originates from the n = term Tq (a, T). This expression 
is different from Eq. (28) in Ref. [2l| . where the Leontovich impedance was used. First, 
the coefficient 1/2(1 + Thomas- Fermi correction) in front of the ^-function shows that only 
p polarization contributes to the A-independent part and, second, the A-dependent part 
behaves as A 2 instead of A In A. These changes are due to different behavior of the exact 
impedance Eq. (jl8j) in comparison with the Leotovich impedance f!24j) . On the contrary, for 
A ^> 1 the Leontovich impedance is a good approximation and we successfully reproduce 
Eq. (33) of Ref. 
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kT 



- a + i- AvpUJa - (Lz^El ] 0( s 

V3cu„ 3^3 V A 



i —2s 



(38) 



where p\ is a numerical coefficient the same as in Ref. [21[, p\ = 0.0133. The only new 
feature in this relation is the presence of the Thomas-Fermi correction. Note that in the 
case of large A both polarization contribute equally to the A independent term (1/2+1/2+ 
Thomas- Fermi correction). 

V. ENTROPY AND DISCUSSION 

Before discussing the entropy behavior in the low temperature limit we should fix the 
parameter a in Eq. (j3J) for the n = term. Let us separate it in two parts describing s and 
p polarizations 

a = a s + a p . (39) 

Contribution of p polarization in the classical part of the free energy JF (a, T) is not problem- 
atic. As we know the only new feature that appeared due to nonlocality is the Thomas- Fermi 
screening. It has clear physical meaning and should present in any reasonable approach. To 
find a p we has to take the impedance (fT9"j) at £ — ► and calculate the function — G p (0)/2 
(see Eq. (J7J)). But we already found the function G p (t), which is given by Eq. (j3*Uj) and in 
our approximation it does not depend on r at all, therefore, 

!^ (4Q) 



2 V V3 c u p 

It is important that the Thomas-Fermi correction in a p is exactly canceled with that in the 
^-independent part of the free energy (j3*7|) or (|3*Hj) . Therefore, the Thomas- Fermi screening 
finally does not contribute to the temperature dependent part of the free energy. 

The real problem is connected with the value of a s . In Bostrom and Sernelius approach 
jlfi| s-polarization does not contribute to the n = term and a s = 0. When SDM prescrip- 
tion is used for the n = term 181 the contribution of s polarization is the same as for 

fl 

the ideal metal: a s = 1/2. The plasma model prescription for the n = term [12[ gives 
a s = a s {uj a /uj p ) as a function of the separation, which approaching 1/2 at u a /u p 1. Close 



13 



value of a s gives extrapolation of the Leontovich impedance from the infrared range to zero 
frequency used in Ref. [25]. Different values of a s are responsible for different temperature 
corrections in these approaches. 

At very low temperature when A ^ 1 the entropy calculated from Eq. (|37|) is 



S 



dT 



k 



?7T(J 



«sC(3) - -0.0146A 2 



(41) 



It goes to zero at T — ■> only if a s = 0. In this case the entropy approaches zero from 



the negative side as T 2 / 3 . This conclusion coincides with that made in Refs. 22 , 23 1 on 
the basis of finite residual resistivity. The use of finite uj t (0) was criticized in Ref. 2f| (see 
a, so a receat preprint Q, on ground that the Ne™. hea t ta was formulated 
for equilibrium states and any defects in the material responsible for the residual resistivity 
should be considered as deviation from equilibrium. The objection is reasonable but here we 
showed that the residual resistivity did not play physical role at low temperatures. Instead 
the nonlocal effects are responsible for the effective relaxation frequency (see Eq. (fTTj) ) 
Vpk ~ Vf/cl ~ 10 13 rad/s, which is much more important than tiny u; r (0). Nevertheless, as 
Eq. (14 lj) demonstrates the final conclusion of Refs. 0, [23 1 holds true. 

At higher temperatures when A ^> 1 but r < 1 the entropy is still negative. In this 
range from Eq. (|3*Hj) one finds for the entropy 



S 



k 



rC(3) 



a, - 



1 64 (l-2p x 

2 + Vl V A 



A> 1, 



(42) 



8vra 2 

which is obviously negative for a s = 0. 

The entropy is a positively defined physical value and the negative value for the Casimir 



entropy is puzzling. Recently [35j some arguments were provided justifying the negative 
Casimir entropy as long as the total entropy is positive. The free energy of the whole system 
consists of two contributions. The main additive part comes from the short-range atomic 
interaction. The long-range interaction realized via fluctuating fields gives much smaller 
contribution to the free energy, but this contribution can be separated due to its nonadditive 
character (see discussion of this problem in Ref. [l^)- The additive and nonadditive parts 
are independent on each other because the first is defined by the volume but the second 
depends on the separation between bodies. Usually it is assumed that the nonadditive part 



is given by the Casimir free energy of fluctuating fields. The idea proposed in Ref. 



is 
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that part of the nonadditive free energy can belong to the bodies. In this case one can write 



AT (a, T) = AT body (a, T) + AT field (a, T) (43) 
Both terms give contribution to the entropy 

J_AT^_d_AT^ d 

v ' ; dT dT K J 

The second term here is negative at low T but the first term could provide the total entropy 
to be positive. This idea can be true but we would like to stress that the term ATb ody (a, T) 
should be explicitly specified. This is because it gives contribution not only to the entropy 
but also to the force according to the relation 

F(a, T) = JA^L _ d A^±. (45) 
da da 

At the moment we do not know any corrections to the Casimir force which appear not from 
the fluctuating field but from the nonadditive free energy of the bodies. In our opinion the 
negative Casimir entropy is the evidence of a thermodynamic problem. We should stress, 
however, that all the other approaches to the temperature correction equally suffer the 
thermodynamic problem because for a s 7^ the entropy is finite at T = 0. 

On the other hand, the zero contribution of s polarization to the n = term has solid 
physical grounds. In the local case the 1/u behavior of the dielectric function, responsible 
for the vanishing of the reflection coefficient r s , is the direct result of the Ohm's law. Any 
attempts to change this behavior will break this law. The plasma model describes well the 
infrared optics but this is only an approximation, which cannot be used as a low frequency 



limit as was proposed in Ref. 



17| . Otherwise any real metal would be a perfect conductor. 



The same is true for the impedance approach, which is extrapolated from the infrared optics 



to zero frequency [25j, [26], and for the SDM prescription as in Ref. [18]. Our nonlocal 
analysis does not bring anything new in the n = term because in the zero frequency 
limit the impedance Eq. ()18j) coincides with the exact local impedance. It is known that 
properly defined local impedances reproduce the force in the dielectric function approach 
3^|, therefore, a s = 0. 

There is a very simple physical explanation why s polarization should not contribute to 
the force in the low frequency limit. If z is the normal direction to the metal surface, then 
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s-polarized field can be chosen as having the following nonzero components of magnetic and 
electric fields: H x , H z , and E y . When u —>■ the magnetic field can be found from the 
Maxwell equation V x H = 47rj/c, where j is the external current density responsible for the 
fluctuating fields The electric field, which is described by the equation V x E = iooH/c, 
will be suppressed in comparison with H because u is small. So in the limit u — > s- 
polarized field degenerates to pure magnetic field. But the magnetic field penetrates freely 
via nonmagnetic metals that means that the reflection coefficient is going to zero. Similarly 
the p-polarized field degenerates to pure electric field in the u — > limit. The electric field 
is screened by the metal and the reflection coefficient is 1. 

We came to a contradictory situation. From electrodynamics it follows that a s = 0. On 
the other hand, thermodynamics shows that the Casimir entropy in this case is negative 
and something must be wrong. All the other approaches proposed in the literature are 
equally unsuccessful thermodynamically (S ^ at T = 0) but, in addition, they do not 
follow from electrodynamics. We cannot resolve the thermodynamical problem by break- 
ing the laws of electrodynamics. Specifically we should stress that the approach based on 
extrapolation of the Leontovich impedance from infrared to zero frequency [3] cannot be 
accepted as physical. It disregards g-dependence of the impedances, which plays crucial role 
for evanescent field configurations. The authors postulated that the evanescent fields have 
the same reflection coefficients as the propagating fields. The Casimir effect is not the only 
physical phenomenon where the evanescent fields can be probed. In the well investigated 
domains like near field optics or near field microwaves q dependence plays principal role. No 
deviations from the standard electrodynamics were noted so far. 

To all appearance the experimental situation is not in favor a s = 0. This case contradicts 
to Lamoreaux experiment |2[. Also there were claims that a s = does not agree with the 
experiments by Decca et al. Q|. However, very high roughness of metallic films in these 
experiments did not allow these claims to be considered seriously. Recently s| the same 
group refined their measurements reducing the surface roughness and increasing precision of 
determination of the absolute separation. It is important that an experimental error of 0.6% 
holds in a wide range of separations from 170 nm to 300 nm. However, in this experiment no 



4l| 



attempt was made to characterize the used gold films optically. Instead, the handbook 
optical data were used for calculation of the force. It was demonstrated that 13 that 
the optical data for gold films prepared in different conditions can variate very significantly. 
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Prediction of the force with the precision better than 2% should include direct measurement 
of the optical properties of the films especially in the mid-infrared range Nevertheless, 
even with the use of the handbook optical data one can conclude that the case a s = 0, 
probably, is not supported by the experiment. This is because the handbook optical data 
present the best samples. The unannealed films used in the experiment should have smaller 
reflection coefficients than the handbook data predict. As the result, the theoretical force 
was overestimated in Ref. JjJ. It means that the difference between the measured force and 
predicted one in the case a s = can be only larger. Of course, this is the result of only one 
group and one has to wait an independent confirmation of it. It should be mentioned also 
that the best way to see the temperature correction is the change of the temperature in 
the experiment. 

All the discussion above shows that the situation with the thermal correction to the 
Casimir force is in deep crisis. To the moment we do not know any approach which is in 
agreement with both electrodynamics and thermodynamics. 



VI. CONCLUSIONS 



We analyzed behavior of the Casimir free energy at low temperatures. The main con- 
tribution to the temperature dependent part of the free energy AjF is defined by the low 
frequencies £ ~ lnkTjfi that is in contrast with the temperature independent part, which 
is defined by the characteristic frequency ( ~ cu a . With the temperature decrease the 
anomalous skin effect becomes increasingly important for AjF. General theory of nonlocal 
impedances was used for calculations. It was demonstrated that at low temperatures the 
relaxation frequency does not play any physical role. Instead, the physical significance get 
the frequency Vpk, where k is the wave number. The approximate Leontovich impedance 
describe the situation well if T ^> (vF/c)(hw a /2ir). When this condition is not satisfied one 
cannot use the approximate Leontovich impedance any more. 

The troubling n = term in the Lifshitz formula was parameterized by the parameter 
a (see Eq. which is different for different approaches to the temperature correction 

discussed in the literature. This parameter was kept arbitrary in calculations. In the tem- 
perature range ftw J v p /c) <C litkT <C hu> p (vF/c) we reproduced for the free energy the same 
result as in Ref. [2l|, where the Leontovich impedance of the anomalous skin effect was 
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used. However, at smaller temperatures, liskT -C hw a (vF / 'c), the behavior of AjF dras- 
tically changes because dependence of the impedance Z s on the transverse momentum q 
becomes important. It was demonstrated that the entropy is going to zero in the limit 
T — > only in the case when s polarization does not contribute to the n = term (a s = 0). 
In all other cases the entropy is finite at T = 0. 

However, even in the case a s = the entropy at low temperatures is negative that, in our 
opinion, indicates the presence of the thermodynamic problem. It was demonstrated that 
the idea on total positive entropy proposed in Ref. (3^ at least incomplete. We concluded 
that the thermal Casimir force is in deep crisis and any approach to resolve the problem 
should respect both the laws of thermodynamics and electrodynamics. 

Note added After submission of this manuscript we became aware of the work by Bo 
Sernelius, Phys. Rev. B (to be published) who also analyzed the nonlocal effects in the 
Casimir problem. The conclusion on the entropy behavior coincides with ours but the 
method of analysis is different. 
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